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Abstract 

Detailed analysis of the distribution of nuclear matter in the halo nucleus 11 Be is performed 
within the microscopic approach based on the flexible variational basis of the functions of the 
"polarized" orbitals model. It is shown that the proposed approach reproduces, in agreement 
with the experiment, the most distinguishing features of the mass distribution in this system. 
Taking into account the results of specific calculations and the exceptionally high reliability of the 
experimental data on the nuclear density distribution in the "tail" of 11 Be, we make a conclusion 
about the existence of new structures with the sizes considerably exceeding the size of the nucleus 
(up to tens of fm) and the density being thousands times less that the typical nuclear one. 

PACS numbers: 21.10.Gv, 21.60.-n, 21.60.Cs, 21.60.Gx 
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I. INTRODUCTION 



As is known, during the last 10-15 years, a number of nuclides with the anomalous N/Z 
ratio has heen discovered (see, for instance, Q, g g Q, 0, Q). These nuclei are 
far from the beta-stability area and may closely approach the boundary of the nuclear or 
nucleon stability. In most of the cases, these are the strongly neutron-excessive isotopes 
of light and middle elements, which include a number of exotic nuclei with interesting and 
unusual properties and structure. Naturally, these nuclides attracted considerable attention 
of both experimentalists and theorists. 

Let us say a few words about the term "halo nuclei". The latter has become very popular 
in recent years and is also used for the well known nuclides (like 6 He, 8 He, etc.). For this 
reason, in order to bring, at least partially, the terminology and the physical meaning into 
agreement, there has been one more extravagant term, "Borromean nuclei", introduced 
Borromean or, alternatively, halo nuclei of Il-nd kind, represent weakly coupled and very 
much fragile nucleon systems being at the same time rather compact in size. They should 
be more exactly referred to as "neutron skin" rather than "halo" nuclei. 

Quite different situation is observed in the case of halo nuclei of I-st kind (the one- 
neutron halo nuclei n Be and 19 C). These nuclides are unusual from every point of view. 
Let us mention the anomalously high cross-sections of the interaction with various nuclei, 
the unexpected inverted spectra of the lower energy levels, and the most unexpected mass 
distribution with the compact core nucleus in the center accompanied by a long sharply 
directed "tail" of very thinned nuclear matter jl|, 

We are not going to consider in detail the basic features of halo nuclei, since there is a 
nu mb e r of exhaustive reviews on the lectin the literature (see, for instance, g g g 
0| , 0). Let us cite only the experimental results on the distribution of the nuclear matter 
in n Be keeping in mind that these latter will be the major subject of our further theoretical 
consideration. These results have been obtained by a large team of Japanese physicists in 
RIKEN |2], where the dependence of the density distribution for the nuclear matter in the 
nucleus n Be on the distance from the nucleus center has been studied. The experimental 
points obtained there have been put within some stripe area bounded by a bold line, in 
order to make allowance for the possible inaccuracies in the determination of density (see 
Fig. 2 in [3]). The authors' very convincing and rather strong statements are cited below. 
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In first place, the existence of the halo tail in n Be nucleus is associated with the last 
weakly coupled neutron only, and, in addition, this valence neutron is the main material 
provider for the "tail" density distribution of nuclear matter. The contributions of the rest 
of nucleons to the "tail" are negligible. The magnitude of the density in the halo tail is 
more than three orders less than the typical density of nuclear matter. 

In second place, while in determination of the central and intermediate (with respect to r) 
regions there is a possibility of effects of additional factors, which, in principle, may bring the 
magnitude of density beyond the area of experimental data obtained by the authors, the slope 
and the magnitude in the "tail" distribution are determined very reliably and accurately. 
The question is, what are the reasons for such confidence of the authors? The point is that 
the most sensitive process with respect to the distribution of the matter in the "tail" is that 
of n Be and 10 Be fragmentation, and the cross section op of this process is anomalously high. 
As a result, the experimental determination of the above mass distribution in n Be can be 
made with such an impressive accuracy. 

One more point important for the theory, which has been emphasized in the paper jsj], 
can be concisely formulated as the indifference of the halo tail in relation to the core nucleus. 
In other words, the presence of the halo tail in n Be almost does not affect the magnitude 
of the basic parameters of the core nucleus 10 Be. 

All the above mentioned relates to the experimental situation on halo nuclei. As regards 
the theory, the mainstream of creative ideas here is focused on the cluster model. There are 
hundreds of papers in the literature, where the successful description of exotic halo nuclei 
is based on the employing one or another of cluster modes, e.g., the core plus one valence 
neutron, or the core plus two valence neutrons. At first sight, such interpretation appears 
very attractive. However, a closer examination reveals here serious difficulties. Let us point 
out some of them. To begin with recall the size parameters of 11 Be. The root-mean-square 
radius of the nuclear matter in n Be is R™ ms = (2, 90 ±0, 05) fm jij, [l^l, whereas the radius 
of the core nucleus 10 Be is R™ ms = (2, 30 ± 0, 02) fm. In order to explain the experimental 
spectrum of the Coulomb decay, the value of the root-mean-square radius of the valence 
nucleon and the diameter of the halo tail must be, in accordance with ~ 7/m, 

and D ha i « 30 fm, respectively. If one assumes that the valence nucleon moves outside 
the core (the concept underlying the cluster representation), a serious question arises as 
to the origin of the long-range interaction between the valence neutron and the nucleons 
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of the core. To treat this problem, one may introduce the so-called "shallow neutron-core 



potential" [12j, and, by appropriately adjusting the parameters for this potential, achieve 
success in calculating specific physical quantities. It is clear that such an interpretation is 
rather phenomenological, and the situation here could be hardly said to be satisfactory. 

The situation with halo nuclei could be clarified within the consistent microscopic 
approach |2] based on the many-particle Schrodinger equation for the wave function 
ty(l,2, A) with some particular choice of the "fundamental" pair NN-potential. Actually, 
the dependence of the wave function \l/(l,2, A) on the coordinates of all A nucleons en- 
ables one to calculate almost arbitrary nuclear property and to reveal the underlying physical 
picture by means of the analysis of relevant theoretical and experimental quantities. 

Taking into account all the above mentioned, we perform in this work exactly this analysis 
for the case of one-neutron halo nucleus n Be. In Section II of the paper, we give a brief 
presentation of the algorithm for calculations within the framework of the polarized orbitals 
model Q| taking an accurate account of the hierarchy of basic principles of the theory. 
In Section III, we present the results of calculations of the nuclear matter density in the 
nucleus n Be. The discussion of theoretical data, along with comparison between the known 
experimental data and similar results of other authors, is given in Section IV. 

II. MANY-PARTICLE WAVE FUNCTION AND THE VARIATIONAL PRINCI- 
PLE 

In a microscopic model of interacting particles, the starting point is the relevant 
many-particle Hamiltonian, where, as a rule, only pair interactions are taken into account, 
i.e. 

A h 2 A 

H A V = E^/; H A = -J2 - T c . m . + £ V itj + U Co uio mb ■ (1) 

Qj=i »>i=i 
, which provides the basis for our further consideration, the pair 

inter-nucleon interaction employed therein is taken in the form of effective central exchange 

NN-potential with the radial dependence specified by the superposition of five Gaussian 

terms [3| 

V « = E E ^ ] + i,2T + i ex P {-^=^£} • PsAhJ), (2) 



where Vjg +1 2T+1 and y, v are the intensity and the radius of the v -th component of the 
interaction potential between the z-th and j-th nucleons, respectively; Ps,t is the known 
projection operator (see, for instance, which separates out the two-nucleon states 

with the spin S and isospin T in the wave function 

Along with the nuclear interaction, the Hamiltonian (1) includes also the Coulomb re- 
pulsion between the protons in the nucleus, which can be conveniently written in the form 
of the integral representation 

Ucoulomb = 2^ |r _ r j = ~J= 2^ / e " (l_ T 2)3/2 ' ^ = ~ Vl ~ ^ 



i>j=l 1 n v i>j=l 



Thus, formally, we deal with the conventional many-particle problem for bound states. In 
view of the possibilities of modern computers, it is worth to seek the solution to the problem 
(l)-(3) on the basis of the well known variational Ritz' principle (see, for instance, [17]). In 
doing so, it is natural to describe the nuclear system in terms of individual particle states 
with the total wave function of the nucleus being approximated by the antisymmetrized 
product of one-particle functions 

V E (1,2,...,A) = -=det[(pi(j)]; <pj(i) = ^(r,-; a h b u d) ■ Xi ms (j>i mt (i), (4) 

where Xi ms and r\\ are, respectively, the spin and isospin wave functions for j-th nucleon 
jis| . and ipn. are the one-particle orbitals. It is clear that the success of such treatment 
will mainly depend on the appropriate choice of variational basis for one-particle orbitals. 
Conventionally, for the functions if) ni , one employs the functions of harmonic oscillator. 
In our model, which is also called in Ref. ^| as "the model of polarized orbitals", each 
single orbital is represented by an oscillator function as well. However, in contrast to other 
authors, we use for every z-th one-particle state its own independent variational parameters 
(di,bi,Ci). 

Note that, as dictinct to the cluster models, all the nucleons equally participate here in 
forming the dynamics of the system. 

The optimal values for the variational parameters {<2j,6j,Cj} as well as the optimal con- 
figuration (determined by the set of quantum numbers {rij = (nf , nf , nf)}) are found as 
a result of minimizing the functional /^({nj; dj, 6j, q}) for the total energy of the nucleus 
under consideration. If the chosen variational basis of the functions ip ni is flexible enough, 
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we can expect that the solution wjj^ a (l, 2, A) to the many-particle problem found in 
such way would be sufficiently realistic. 

We are not going to dwell upon the so-called "technical" issues, although, exactly these 
latter frequently determine the success of a work. Let us only list the steps made for eval- 
uating the functional /^({n^; a«, bi, Cj}) :i) matrix elements for the operators of all physical 
quantities entering the Hamiltonian (1) have been obtained; ii) projection of the internal 
wave function (4) onto the energy state J n = \ + , associated with the ground state of 11 Be 
has been carried out on the basis of the Hill- Wheeler integral; iii) search for the global mini- 
mum of the functional /^({nj; cij, bi, q}) for the total energy of the nucleus in the hyperspace 
of variational parameters {at, bi, c{\ has been performed. 



III. RESULTS OF NUMERICAL COMPUTATIONS 

As mentioned above, the calculations of the needed wave function a (l,2, .... A) 
and some basic features of the nucleus n Be have been performed in previous works [14]. 
Therefore, we can use here the results of jl^ in order to evaluate the sought-for distribution 
of the nuclear matter in n Be. 

Recall that, according to [14], the ground state of the nucleus n Be is associated with the 
configuration: 

protons Z = 4 -> (000) 2 (001) 2 ; neutrons N = 7 -> (000) 2 (001) 2 (010) 2 (002) 1 , (5) 

i.e., the protons have the regular subsequently filled configuration, whereas the neutron 
configuration is unusual. The optimal values for variational parameters of the orbitals il> ni 
are given in Table 1. 

Note that the physical meaning of the variational parameters (aj,6j,Cj) is that they de- 
termine the extent of spreading for the i-th one-particle state along x, y, z axis, respectively. 
The neutron configuration indicated in (5) turned out to be energetically favorable only in 
our model of "polarized" orbitals, where for the first time in calculations of such kind the 
variational parameters {dj, bi, c{\ were independent for any of the orbitals ip ni and in any of 
the directions x, y, z. Unproportionally large value of the variational parameter c!\ = 5.38 fm 
(see Table 1) evidences the anomalous spreading of the valence neutron along one of the 
coordinate axis. 
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The information on the explicit expression for the wave function enables us to evaluate 
the sought-for mass distribution qa{t) in the nucleus 11 Be. To this end, we have to average 
the operator for the one-particle spatial density with the wave functions \l/^- a (l, 2, A) , 
i.e. 

A 

QA(r) = (^ I jY, 6 ( r ~ r ^L)- ( 6 ) 

i=l 

Here the coordinates are written in center mass system R c . m . = 0. 

Note that in experiments, only the radial dependence for the nuclear density distribution 
averaged over the angles fl 

I p2tt pit 

Qa{ t ) = k / dd I qa(t) sm8d8 ; r > Ifm (7) 

47rr 2 J Jo 

is measured. In the vicinity of the point r = 0, the averaging over the angles Q has to be 
performed as 



N, 



point 



1 

Q A {r) « — ; r < 2 ^ m ' ( 8 ) 

^points 

where the points are taken on the sphere of the radius r = |rj|. 

The distribution of nuclear matter for the wave function of the ground | + state of the 
nucleus calculated in the above way is shown in Fig.l. Actually, this figure represents an 
almost identical copy of Fig. 2 of Ref. Q], where, in addition, our theoretical curve (the line 
with dots) is displayed. The three curves with /? = 0, j3 = 0.5 and (3 = 0.7 in Fig.l are the 
attempts of other authors to describe the experimental situation within the framework of 
the deformed shell model (see the references in 0). As is seen from Fig.l, the theoretical 
data (the line with dots) and the experimental estimates for the nuclear density distribution 
in the nucleus 11 Be are in good agreement. It is particularly true for the magnitude and 
the slope in the "tail", where the experimental estimates are very reliable and a complete 
coincidence between the theory and experiment is observed. 

It is worth mentioning that in theory, in contrast to experimental studies, there exists 
a possibility to investigate the above mass distribution in the internal coordinate system 
rigidly coupled to the nucleus. In this case we deal with the spatial structure Qa{t), which 
is the function of three coordinates x, y, z. In Fig. 2, we give one of the sections (in (x, z) 
plane) for the function qa{y) for the nucleus 11 Be under consideration. The figure in the 
plane (y, z) would be similar, as far as the size properties are concerned. The choice of the 
values for the lines of constant density in Fig. 2 has been made in agreement with the relevant 



experimental data within the interval [5 fin < r < 12/m]. Finally, in Fig. 3 we present the 
distribution of nuclear matter Q neu tron{ r ) for the valence nucleon, or, more exactly, the more 
exactly, the lines of constant density in the plane (x,z). In this figure, we clearly observe 
three structures, the two lateral ones, symmetrically placed at the distance of 10/m, and the 
central one being considerably smaller in size than the lateral structures. All these three 
structures have approximately equal density with the magnitude (2 — 5) ■ 10 _5 /m -3 . 

IV. ANALYSIS AND CONCLUSIONS 

After introduction to both experimental and theoretical results on the mass distribution 
in halo nucleus n Be, let us try to make analysis in order to obtain a physical interpreta- 
tion. Let us start with Fig.l. The radial distribution for the density of nuclear matter 
presented here could be hardly brought into agreement with the notion of "atomic nucleus" 
conventionally given in reference books. Specifically, the nucleon density is said therein to 
be nearly constant in the nucleus center, and to decrease exponentially at periphery. As is 
seen from Fig.l, the nucleus n Be is not the case. Other features of nucleon system under 
consideration are in contrast with the classical concept of "atomic nucleus" as well, let alone 
the size properties. 

Fig. 2 and Fig. 3 strongly evidence that in the case of the halo nucleus n Be we deal, in 
fact, with two distinctly different physical systems, the well known nucleus 10 Be referred here 
to as the "core nucleus", and the structure of intermediate, between atomic and nuclear, 
nature. The latter is presented by three separate objects. This structure is, first of all, a 
result of the action of the Pauli principle with respect to all 11 nucleons constituting the 
11 Be system. It is coupled to the core nucleus by its central part, which completely envelopes 
the core nucleus. Two lateral much more massive parts are outside the nucleus 10 Be, i.e., 
beyond the action of nuclear forces. These three parts-objects are held together only due to 
the fact that they originate from a single quantum of nuclear matter, the neutron. 

To put it in the other way, the nuclear processes give rise to the conditions, which make 
it energetically favorable for the system of 4 protons and 7 neutrons to form the compact 
nucleus 10 Be and an unusual structure with the mass and the quantum numbers being 
exactly equal to the valence neutron. According to the experimental estimates, this energy 
gain is equal to only 5E = 0.32 MeV. Exactly at this distance, on energy scale, the first 
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excited \ state of the nucleoli system under consideration is observed. 

Thus, the analysis of the experimental situation on the distribution of nuclear matter 
in the exotic nucleus n Be performed in this work has shown that the nucleon system n Be 
should be considered as two different in nature subsystems, the compact nucleus 10 Be, and 
some structure emerging in the vicinity of this nucleus as a result of Pauli principle, which 
is distinctly different in its features from an atomic nucleus. Therefore, we should regard 
the discovery of the halo nuclei of I-st kind as the discovery of new structures of microcosm 
with the sizes ranging up to tens of fm and the density being thousand times less than the 
nuclear one. 
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TABLE I: The optimal values of the variational parameters for one-particle orbitals ifj ni 



Parameters {aj,6j,Cj}. 


The optimal values, fm. 


V 

°i 


1.355 


V 
«2 


1.49 


6? 


1.55 


fig 


1.685 


c? 


1.84 


eg 


1.92 


a? 


1.375 


4 


1.51 




1.61 


a n 


1.725 


b n 


1.52 


°2 


1.64 


°3 


1.80 


h n 
°4 


1.90 


r n 
C l 


1.89 


r n 
c 2 


1.96 


c 3 


2.18 


r n 

c 4 


5.38 
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FIG. 1: Density distribution of nucleons on logarithmic scale for the ground ^ state of Be. The 
results of RIKEN : 8] are shown by thick lines; the results of our calculation are shown by thin line 
with filled circles. The thin curves with f3 = 0, (3 = 0.5 and (3 = 0.7 are the results of the deformed 
model (see Jsj]). 
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FIG. 2: Intrinsic structure of ^ state of Be. The figure shows sections of the matter density 
qa(x,u = 0,z) in the plane (x, z). 
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FIG. 3: The density distribution for the valence neutron in the plane (x, z). 
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